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Abstract 

There are many materials whose dielectric properties are described by a stretched 
exponential, the so-called Kohlrausch- Williams- Watts (KWW) relaxation function. 
Its physical origin and statistical-mechanical foundation have been a matter of de- 
bate in the literature. In this paper we suggest a model of dielectric relaxation, 
which naturally leads to a stretched exponential decay function. Some essential 
characteristics of the underlying charge conduction mechanisms are considered. A 
kinetic description of the relaxation and charge transport processes is proposed in 
terms of equations with time- fractional derivatives. 
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1 Introduction 



A rich variety of materials with molecular or structural disorder have two 
important properties in common. One is non-exponential, non-Debye character 
of dielectric relaxation, often described by a stretched exponential, the so- 
called Kohlrausch- Williams- Watts (KWW) function [H2] 



Mt) = eM~(t/rf] (1) 
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with the exponent < (3 ^ 1 and r a constant. This stretched exponential 
relaxation has been found empirically in various amorphous materials as for 
instance in many polymers and glass-like materials near the glass transition 
temperature (for reviews see Refs. [3] and [I], and references therein). Some 
physical models incorporating general features of the stretched exponential 
dielectric relaxation are discussed in Refs. [HI51I6] where one also finds reviews 
of experimental dielectric relaxation data. 

The other property, often found in disordered insulators and semiconductors, 
including those doped with electric charge, is universality of ac (alternating- 
current) conduction, which is expressible in terms of a power-law dependence 
of the real part of conductivity on frequency (Ref. [7] for review): 



with the exponent ^ r\ < 1. The universality means that the r\ values 
do not depend on the details of the underlying conducting lattice nor on the 
microscopic charge transport mechanism operating in the system (i.e., classical 
barrier crossing for ions and/or quantum mechanical tunneling for electrons). 
Both (3 and 77 depend on the chemical composition of the material and the 
absolute temperature. 

The physical background and the statistical mechanical foundation of the 
above originally phenomenological expressions have been a matter of debate 
in the literature. Here we mention two hard results in their favor: (i) The 
KWW function expands into a weighted superposition of the Debye single 
exponential relaxation functions |5J 



where the weighting function Qp{[x) is expressible in terms of a stable (Levy) 
distribution. Due to this connection with the statistics of stable laws, stretched 
exponential relaxation functions can be argued to appear naturally, thus be- 
ing a characteristic property of systems in which the dynamics occur on many 
time scales, (ii) The power-law form of ac conduction coefficient derives from 
a model [Sj, in which the conduction occurs as a result of random walks of 
charged particles on a percolating cluster. In this paradigm the universality 
of ac conduction is rooted in the universality of the percolation transition [9]. 
Support for this standpoint can be found in the results of Ref. [10] . The expo- 
nent T) is expressible in terms of the percolation indices [8] or the topological 
characteristics of the lattice [H] . 
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In this study we analyze the properties of dielectric relaxation in disordered 
solids on an equal footing with the ac conduction properties. In section 2 we 
demonstrate that the exponents (5 and rj are related to each other via (5 = l—rj. 
We then suggest simple numerical estimates of the r\ values in a regime when 
the conduction concentrates on a percolating cluster. In section 3 we formu- 
late a self-consistent model of dielectric relaxation, in which the KWW decay 
function is obtained from a power-law dependence of ac conduction coefficient 
on frequency and the basic electrostatic equations. In section 4 we describe 
the KWW relaxations kinetically. We are led to the issue of kinetic equations 
with time fractional derivatives and we propose a systematic derivation of the 
fractional relaxation and fractional diffusion equations from the property of 
ac universality. We summarize our results in section 5. 



2 Stretched-exponential relaxation and ac universality 

2. 1 General 

Let a homogeneous, isotropic dielectric be exposed to the external polarizing 
electric field E = E(t, r), which we consider as a function of time t and the 
space coordinate, r. By homogeneous and isotropic we refer to spatial scales 
larger than the typical scales of the molecular /structural disorder. Assuming 
a linear and spatially local response of the material the polarization field at 
time t at point r can be written as 



where x{t ~ t') is a response or memory function. Causality requires that 
x(t — t') = for t < t'. A Fourier transformed memory function x{t) is defined 
to be the frequency- dependent complex susceptibility of the material, i.e., 



Given x(^) one defines the frequency-dependent complex dielectric parameter 

as 




(4) 



— oo 




(5) 



— oo 



e(u) = 1 + 4ttxH 



(6) 
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In a basic theory of the decay of polarization one is interested in the dielectric 
response to a field which is steady for t < and, then, is suddenly removed 
at time t = 0. It is then common to write the complex function e(u) as [TTJ 



e(w)-e(oo) ld(j>{t)^ 



e(0) - e(oo) J dt 



dt (7) 



where e(oo) and e(0) are, respectively, the high and zero frequency limits of 
e{u) and (f)(t) is the function that describes the decay of polarization after 
the polarizing field has been removed. Equation (J7|) is essentially an integral 
equation for the decay function. The decay function can be obtained directly 
from the magnitude of the polarization response to an electric field with time 
history as discussed in Ref. [12]. 



2.2 Connection between relaxation and memory functions 



It is easy to obtain the connection between the functions <p(t) and x{t)- We 
have, with = H^) — 1] / 47r, 



1 e(0)-l ^U) 1 



X (t) = — / -H e~ wt du 

AK ' 2tt 4tt J e(0) - 1 



-oo 



In writing Eq. (jSJ) we took into account that and x(°°) form a Fourier 
pair. Assuming that, in the high frequency limit, the value of the dielectric 
parameter is, approximately, e(oo) = 1, and combining Eqs. ([7]) and (jHJ) we 
find 



where i = e(0) — 1 is the dielectric loss strength. Changing order of integration 
in Eq. ([9]) we write 

-oo 



Considering that /+~ e'^-^dcu = 2nS(t - t') from Eq. QIQ} it is found that 
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Hence the memory function x{t) is just the time derivative of the relaxation 
function <p(t) (with a proper coefficient in front). Setting <fi(t) = exp[— (£/r)^] 
in Eq. (fTTj) we get 

x{t)=(3^-/- 1 eM-(tlrf] (12) 



Note that x(t) diverges as W 1 for t — ► +0 when (3 < 1. 
2.3 Leading-term approximation 



We now turn to an explicit calculation of the o~'(u) function from a KWW 
decay function. In fact, setting </>(£) = exp[— (t/r) 13 } and integrating by parts 
in Eq. (J7J) we find, for the real and imaginary parts of the complex dielectric 
parameter, 

e'(u) - 1 = £[1 -ttzV(z)] (13) 
e"(u) = hzQ(z) (14) 

where z = ujt is the dimensionless frequency, u = t/r the dimensionless time, 
and V(z) and Q(z) are the standard integrals: 

1 °° 

Viz) = — I exp( — it^) sin(zu)du (15) 

7T J 

o 

1 °° 

Q(z) = — exp( — uP) cos(zu)du (16) 

7T J 



which play a prominent role in the theory of stable distributions [5]. In the 
higher frequency range, physically corresponding to a departure from the "in 
phase" conductivity response [10] and the transition to ac charge transport 
mechanisms, series expansions of these integrals become [3T5TT3] 

ZV{Z) = n ^T^TlT C ° S ~2~ 

V ; vr;^ v ; r(n + l) 2 v ; 



From Eqs. (ITT)) and ([TBI) one can see that the expansion of ttzQ(z) starts from 
a term which is proportional to z~^ and so does the expansion of 1 — nzV(z). 
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Writing the susceptibility function as 

*H = ^[1 - Tr^(^) + zttzQ(z)] (19) 



then going for the leading term in the power expansion oil — ^zV{z) + mzQ{z) 
we find 

XH oc w-" (20) 



Hence 

eV) - law^ (21) 
e"(w) oc uT^ (22) 

Noticing that the real part of the frequency-dependent complex conductivity 
cr(u;) is expressible as 

a\uj) = uje"(uj)/Atx (23) 

from Eqs. flU]) and (J23]) we get 

= ^Q(z) (24) 

with the leading term 

a'{u) oc cu 1 ^ (25) 

Equation (j25|) reproduces the phenomenological expression in Eq. (j2J) with the 
power exponent 

71 = 1-/3 (26) 



If the conduction occurs on a percolating cluster, the rj values can be obtained 
from [9] 

, = ^ (27) 



Here, # is the exponent of anomalous diffusion |8j, which also appears in de- 
scribing the dc conductivity near percolation, d is the topological (integer) 
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dimension of the ambient real space, and dj ^ d is the Hausdorff fractal di- 
mension of the subset, on which the conduction is concentrated. Using known 
estimates [T4] for 9 and df one has 77 ~ 0, 0.3, and 0.6 for d = 1, 2, and 3, 
respectively. The mean-field result, holding for d > 6, is r\ = 1. From Eqs. (1261) 
and (J2"T|) it is found that 

■ 8 = ^TT^ (28) 



yielding /3 ~ 1, 0.7, and 0.4 for d — 1, 2, and 3. Interestingly, the estimate 
/3 ~ 0.4, holding for d = 3, fits the original result due to Kohlrausch (Refs. 
PP and [1] for review). This value corresponds to the value rj = 1 — (3 ~ 0.6, 
which nicely falls into the range of theoretical arguments [9] and the existing 
observational evidence I71I10II15I. 



3 Self-consistent dynamic-relaxation model 



3. 1 Formulation of the model 



Our goal now is to obtain the KWW decay function analytically. We start 
with a polarization field of the form 

+00 

P(t, r) = P(0, r) + J X (t ~ t')E(t', r)dt> (29) 


where P(0,r) is the initial polarization, and E(t',r) is the electric field. We 
shall assume that there are no external electric fields acting in the system after 
time t = 0, so that E(i', r) is essentially the inherent, self-consistent electric 
field in the bulk of the dielectric due to the electric charges present. 

Let p(t, r) be the density of the electric charges at time t at point r. The 
function p(t, r) is defined as the mean density of the charges in a small volume 
around r such that the highly fluctuating molecular densities are averaged out. 
If there are no external charges, the density p(t, r) is essentially the density 
of the polarization charges, which we shall denote by p(t,r). If the external 
charges are present in the bulk of the dielectric with the density p ex t{t, r) then 

p(t,r) = p{t,r) + p cxt (t,r) (30) 



By external charges we mean charges of external origin, that are alien to the 
material. These charges may be present due to a charge injection. By allowing 
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for external charges we include a class of doped insulators and semiconductors, 
which offer a challenging set of fundamental problems, as well as prominent 
technological applications (Refs. [TBlf!?] and references therein). Equation (IHUi) 
extends the model in Ref. [12] in which p ext (t, r) = 0. 

In the basic theory of dielectrics one writes 

V • E(t, r) = 4vrp(t, r) (31) 

and 

V-P(t,r) = -p(t,r) (32) 
so that 

V-D(t,r) = 4vrp ext (t,r) (33) 

where D = E + AnP is the electric displacement in the medium. The density 
of the polarization currents is defined by 

j(t,r) = ^P(t,r) (34) 
Performing d/dt on Eq. ( 1291) it is found that 

+oo 

j(t,r)= J a(t-t')E(t',r)dt' (35) 
o 



where we introduced 

a(t-t>) = ^ x (t-t>) (36) 



Note that, due to causality, a(t — if) =0 for t < t' '. 

If external charges are present in the system, they may cause, in addition, their 
own current, j ex t(^ r )- Before we write an expression for j ex t(^ r ) we address 
the issue of the microscopic charge transport mechanism: 

It is generally believed that, in relatively poor conductors such as major glasses 
and polymers, the polarization current is caused by orientational motion of 
polar molecules or dipoles containing parts of these [3~][Tl~] . The external cur- 
rent, in its turn, could be thought of as arising from a migration of charged 
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particles along the underlying stationary molecular distribution. This migra- 
tion could be mediated by a specific chemical composition of the material, as 
for instance by the bonding structure in conducting polymers [H]. In porous, 
nano-crystalline materials, the transport of charge may also show a strong, 
nonlinear dependence on the charge density and injection, a phenomenon usu- 
ally explained in terms of trap-filling (Ref. [17] and references therein). As a 
model approximation, here we shall rely on the hypothesis of trap-controlled 
conduction and diffusion, in which the transport occurs as a result of hopping 
[7j of charged particles between the localized states. If the hopping has a char- 
acteristic time, then the transport is described by a Markovian chain process 
with a characteristic hopping frequency. In a more general situation there is 
a distribution of waiting or residence times between the consecutive steps of 
the motion and the Markovian property is invalidated. The current density is 
then a flow with memory: 

+00 

j ext (t, r) = J a cxt {t - t')E(t', v)dt' (37) 

Q 

where (T ext {t — f) is a memory function which describes the multi-scale trap- 
ping and detrapping of the external charges in wide-gap potential wells of the 
conduction-band level. Due to causality, cr ex t(^ — t') = for t < t'. 

The total current density in the bulk of the material can now be written as 

j(t,r)=j(t,r)+j ext (t,r) (38) 

Utilizing Eqs. (1351) and fl3T|) we have 

+00 

j(t,r)= J a(t-t')E(t',r)dt' (39) 




with 



a(t-t') = a(t - t 1 ) + a ext (t - t') 



(40) 



A Fourier transformed a(t) is defined to be the frequency-dependent complex 
conductivity of the material, i.e., 

a(u) = a(u) + a cxt (uj) (41) 



where we introduced the partial ac conductivities ct{uj) and a ext (u) due to 
respectively the polarization and external charges. The conservation of the 



9 



electric charge is expressed by the continuity equation 

^p(t,r) + V-j(t,r)=0 (42) 



Substituting j(t, r) from Eq. (I39I) and taking V- under the time integration we 
have 



0_ 

dt 



+00 

p(t,r)+ J a(t -t')V • E(t',r)dt' = (43) 



Utilizing Eq. fl3Tl) we finally arrive at a closed integro-differential equation for 
the charge density, i.e., 



o +00 
a 



-p{t,r)+4n J cx{t - t')p{t', r)dt' = (44) 



3.2 Stretched exponential relaxation functions 



By Laplace transforming Eq. (I44j) we find 

sp(s, r) - p(0, r) + 47ra(s)p(s, r) = (45) 

where p(0, r) is the density of the charges at time t — 0, and cr(s) is the Laplace 
transform of o~(t). We now speculate on the form of the o~(s) function: 

There is an increasing belief [18] that scale-invariance and fractality are hall- 
marks of chaos and disorder. That such an argument leads to a power-law 
behavior of ac conduction coefficient was pointed out by Milovanov and Ras- 
mussen [5] who based this on the early work of Gefen et al. [8J. One would 
expect that, for those time scales on which the dynamics are dictated by the 
disorder, the partial ac conductivities can be modeled by power-laws, i.e., 
oc u Vl and cr ext (cG>) oc oj m with some fractional r\\ and r) 2 . Here we assume 
that r)i ~ i] 2 , i.e., that the two exponents are approximately the same. This 
assumption refers to universality of ac conduction [7]fTU] and indicates that the 
scaling properties of conduction of polarization charges and external charges 
are determined by the fractal geometric properties of the material and not by 
the details of the conduction mechanism. Writing a{uS) oc and <7 ext (u;) oc u> v 
with the same power rj we have, for the total conductivity 

o-(s) = as v (46) 
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with a a constant coefficient. This power-law form is just the Laplace version 
of Eq. (|2j). Note that a may generally depend on the average concentration of 
the external charges in the conducting domain. 



Separating variables in Eq. we write p(s, r) = <f>(s)ip(r) with ip(r) = p(0, r) 
the initial charge- density. Combining Eqs. (|45p and (|46p it is found that 



1 



s + Arras 11 



(47) 



In the time domain, 



i 



,,Sf 



27ri J s + r p s 

— ioo 



ds 



(4* 



where we introduced the notations /3 = 1 — rj and r _/3 = Ana. Equation 
coincides with the definition of the Mittag-Leffier function E^ —(t/r) 13 (Eq. 
(B.l) in Appendix B of Ref. [IS])- The Mittag-Leffier function has the series 
expansion 



E, 



] r(n/3 + l 



(49) 



For short times, this expansion goes as a stretched exponential, 



<t/r) 



exp 



{t/rf/r(f3 + 1) 



(50) 



This closed analytical form replicates the KWW decay function in Eq. (CQ). 
Finally, for the charge relaxation by ac charge transport mechanisms, 



p(t,r)ocexp \-(t/T)P/r(J3 + l) 



(51) 



where we omitted the space dependence for simplicity. Limiting cases of ex- 
pression ( 1511 are the following: 

p(t, r) ^> pc X t(t, r) — Equation (loTT) leads to a stretched exponential relaxation 
of the polarization charges. For many years, theoretical justification of this 
stretched exponential relaxation regime has been an issue in the theory of 
dielectric relaxation |20j . 

p ext (t, r) ^> p(t, r) — Equation (!5~T|) leads to a stretched exponential relaxation 
of the external charges. This regime may be appropriate for poorly polarizable 
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disordered media when the above condition is satisfied due to a small p(t, r). 
An observational verification of this regime might constitute an experimental 
challenge. 



4 Fractional kinetic equations 



The purpose of this section is to describe the KWW relaxations kinetically. 
We intend to demonstrate that the power-law dependence of the ac conduc- 
tion coefficient on frequency leads to a fractional extension of the relaxation 
and diffusion equations, which accommodate fractional-order time derivatives. 
These offer a suitable analytic formalism to incorporate the features of ac uni- 
versality and the underlying structural disorder. The discussion below draws 
on the paradigm of fractional kinetics |21j . which finds expanding applications 
in various fields of research 22.23,24.251. 



4-1 Fractional relaxation equation 



First, we notice that the power-law s v with the fractional < rj < 1 is 
the Laplace transform of the Riemann-Liouville derivative, which is defined 
through [27] 

1 3 



° D ^ = W^w t J dt 'f^f, (52) 







with ip(t, r) a function from the class of differintegrable functions. The Riemann- 
Liouville derivative is a well-defined fractional extension of the ordinary partial 
time derivative. Setting er(s) = as'' in the dispersion relation (1451) and replac- 
ing s^ by oD? we write, with r\ = 1 — (3 and r _/3 = 47ra, 

^-p(t, r) = -T-P oD^t, r) (53) 



Equation flo^l) is the canonical form of the fractional relaxation equation |28j . 
Applications of this are reviewed in Refs. [T9|23H24j . Here we add to the existing 
knowledge by proposing that relaxations in disordered solids are described 
by the fractional relaxation equation, provided that the dynamics are self- 
consistent, and the property of ac universality is verified. 
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4-2 Fractional diffusion equation describing sub-diffusion 



The above analysis applies to length scales much longer than the mean-free 
paths of charges participating in the ac conduction processes. At length scales 
comparable to or shorter than these, the dynamics of relaxation should be 
described kinetically. The key issue is the form of the flow function, which we 
define as 

t 

i(t,r) = -Jv(t-t')Vp(t',r)dt' (54) 
o 



with a memory kernel Dit — t') such that T>(t — t') =0 for t < t' . When 
V(t — t') is a delta function, Eq. (1541) reduces to the well known, Fick's law. 
A Fourier transformed T>(t) is defined as the frequency-dependent complex 
diffusion coefficient, V{oj). The value of V{uj) can be expressed in terms of the 
ac conduction coefficient as 



T 



ne A 



(55) 



Here, e denotes the carrier charge, n their number density, and T the absolute 
temperature. Equation ( l55l) indicates that the properties of charge conduction 
and diffusion are determined, for each time and frequency scale, by essentially 
the same collisional properties. For the diffusion on fractals, Eq. ( 1551) can be 
obtained as a Fourier transform of the average size- and time-scale depen- 
dent diffusion coefficient [8]. In the zero-frequency limit, Eq. (1551) reduces to 
the conventional Einstein relation between the diffusion constant and the dc 
conductivity. 

Combining Eqs. ([2]) and (1551) we can propose that, in the frequency range in 
which the ac conduction coefficient can be modeled by a power law, 

V(u) oc oo v (56) 



Such power-law behavior of the frequency- dependent diffusion coefficient has 
been found in, for instance, stochastic Hamiltonian systems (Ref. |26j and 
references therein). 

By Laplace transforming Eq. ( 1541) we get 

j(s,r) = -P(s)Vp(s,r) (57) 
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where T>(s) is the Laplace transform of T>(t). When substituted into the con- 
tinuity Eq. (|4"2l this yields 



sp(s,r)-p(0,r)=P(s)V 2 p(s,r) 



(5f 



Adhering to the power-law form cr(s) = as v from Eq. (I55I) we have P(s) = As 77 
with A = aT/ne 2 . Utilizing the scaling P(s) oc s 7 ' in Eq. (1581) we write 



sp(s, r) — p(0, r) = s v V 2 p(s, r 



(59) 



where A = 1 for simplicity. In the time domain, Eq. (|59~I) becomes 
5 



p(t,r)= D t 1 - /J V 2 p(t,r) 



at 



(60) 



where we used r\ = 1 — /3. Equation ( 1601) is the canonical form of the fractional 
diffusion equation describing sub-diffusion [19], with f3 the fractal dimension 
in time [29]. In various settings, this equation has been derived and discussed 
in the literature [TM9T23ll24T25^ 



The characteristic function which is the two sided Fourier transform of p(t, r) 
over the space variable r satisfies 



0_ 

Of 



p(t,k) = -k\Dl^p(t,k) 



(61) 



Equation (16T|) is essentially the fractional relaxation equation in wave-vector 
space. For the small Ak 2 t^ < 1, the characteristic function reduces to a 
stretched exponential decay function: 



p(t,k) ?»exp -kV/r(/3 + l) 



(62) 



In the real space, the fundamental solution of the fractional diffusion Eq. (160!) 
is expressible in terms of a stretched Gaussian distribution (see Ref. [21] where 
further particularities of the initial conditions for time-fractional equations are 
discussed) . 



5 Summary 



We have discussed the properties of dielectric relaxation and ac (alternating- 
current) conduction in disordered solids, treating them on essentially the 
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same footing. Having assumed the property of ac universality, we found that 
the relaxations are stretched exponential rather than the Debye exponential. 
Our results comply with the classical phenomenological expressions due to 
Kohlrausch, Williams, and Watts (KWW). We have shown that the KWW 
decay function can be obtained analytically from a self-consistent model of 
dielectric relaxation, in which both the polarization and electric source fields 
are self-consistently generated by the residual charge-density. The exponent 
of the KWW decay function is related to the exponent of the ac conduction 
coefficient via /3 = 1 — 77. Assuming that the conduction concentrates on a per- 
colating cluster we found r\ values within the range of observational evidence. 
Finally, we found that the relaxations are described by a fractional exten- 
sion of the relaxation and diffusion equations, which naturally incorporate the 
power-law dependence of ac conduction coefficient on frequency. 
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